arXiv:1504.07639v2 [astro-ph.GA] 30Jul2015 


Mon. Not. R. Astron. Soc. 000, 1-15 (2015) 


Printed 31 July 2015 (MN WT^X style file v2.2) 


The First Population II Stars Formed in Externally Enriched 
Mini-halos 

Britton D. Smith^*, John H. Wise^, Brian W. O’Shea^f, Michael L. Norman^, 
and Sadegh Khochfar^ 

^ Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK 

^Center for Relativistic Astrophysics, Georgia Institute of Technology, 837 State Street Atlanta, GA 30332, USA 

^Department of Physics & Astronomy, Michigan State University, East Lansing, MI 48824, USA 

"^Center for Astrophysics and Space Sciences, University of California at San Diego, La Jolla, CA 92093, USA 


31 July 2015 


ABSTRACT 

We present a simulation of the formation of the earliest Population II stars, starting from 
cosmological initial conditions and ending when metals created in the first supemovae are in¬ 
corporated into a collapsing gas-cloud. This occurs after a supernova blast-wave collides with 
a nearby mini-halo, inducing further turbulence that efficiently mixes metals into the dense 
gas in the center of the halo. The gas that first collapses has been enriched to a metallicity of 
Z^2 x10“^Zq. Due to the extremely low metallicity, collapse proceeds similarly to metal- 
free gas until dust cooling becomes efficient at high densities, causing the cloud to fragment 
into a large number of low mass objects. This external enrichment mechanism provides a plau¬ 
sible origin for the most metal-poor stars observed, such as SMSS J031300.36-670839.3, that 
appear to have formed out of gas enriched by a single supernova. This mechanism operates on 
shorter timescales than the time for low-mass mini-halos (M ^ 5 x 10^ Mq) to recover their 
gas after experiencing a supernova. As such, metal-enriched stars will likely form first via this 
channel if the conditions are right for it to occur. We identify a number of other externally 
enriched halos that may form stars in this manner. These halos have metallicities as high as 
0.01 Zo, suggesting that some members of the first generation of metal-enriched stars may be 
hiding in plain sight in current stellar surveys. 

Key words: cosmology, hydrodynamics, radiative transfer, methods: numerical, galaxies: 
star formation 


1 INTRODUCTION 

The physical processes relevant to star formation have conspired to 
produce stars of predominantly low mass for most of the history of 
the Universe. The stellar initial mass function (IMF) appears to be 
so robust to variations in environment that debate over its univer¬ 
sality primarily concerns whether the statement “the IMF is uni¬ 
versal” should include the word almost. Excellent reviews of this 
subject have been provided by Kroupa (2002) and more recently, 
Bastian et al. (2010). The one notable exception to the universal¬ 
ity of the IMF is the case of metal-free (Population III or Pop III) 
stars. The mid 2000s was a boom for quality reviews of this subject 
(Barkana & Loeb 2001; Bromm & Larson 2004; Ciardi & Ferrara 
2005; Glover 2005; Ripamonti & Abel 2006; Bromm et al. 2009). 
In short, neutral metal-free gas cools very inefficiently as it col- 
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lapses, resulting in Jeans mass-scale fragments that are of the order 
of thousands of M©. Early simulations found that this resulted in 
the formation of only one (Abel et al. 2002; Bromm et al. 2002; 
Yoshida et al. 2006; O’Shea & Norman 2007) or two (Turk et al. 
2009; Stacy et al. 2010) dense cores that were free to accrete from 
this massive envelope and grow uncontested. More recent simu¬ 
lations have found that the disks surrounding the massive, central 
core can be unstable to fragmentation, potentially providing a chan¬ 
nel for lower mass Pop III stars (Clark et al. 2011 ; Greif et al. 2011 ; 
Stacy & Bromm 2014). However, simulations that follow the late¬ 
time accretion onto the central object find that they will grow to 
be 10s to 1000s of M© in final size (Hosokawa et al. 2011; Hirano 
et al. 2014, 2015; Susa et al. 2014), making it quite clear that the 
Pop III star formation mode is, at a bare minimum, unquestionably 
distinct. Even on the low mass end of the Pop III IMP, Hartwig 
et al. (2015) claim that the existing sample size of low-metallicity 
star surveys in the Milky Way can already rule out the existence of 
metal-free stars below 0.65 M© with 95% certainty, barring enrich¬ 
ment via accretion of metals from the interstellar medium (Johnson 
2014). 
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Following naturally from the idea that the first stars to ever 
form in the Universe represent an exception to the universal stel¬ 
lar IMF is the question of how this necessary transition in star 
formation modes took place. The Universe was forever changed 
when supernovae from the first stars created the first metals, whose 
introduction to star-forming gas enhanced its ability to cool and 
fragment. From this argument of increased cooling and fragmenta¬ 
tion, two “critical metallicities” have been found for yielding a sig¬ 
nificant departure from the thermal evolution of collapsing metal- 
free gas. At low densities^, ~ cm“^, fine-structure lines of 

atomic C and O allow collapsing gas to continue to cool past the 
limitations of metal-free coolants (primarily H 2 ) at a metallicity of 
~ 10-3.5 (Bromm & Loeb 2003). At much higher densities, 

~ IO 12-14 Oust emission can induce a sharp cooling phase 

when the metallicity is only ~ 10“^'^ Z©, assuming the dust- 
to-gas ratio scales with the gas-phase metallicity (Omukai 2000; 
Omukai et al. 2005; Schneider et al. 2006). Three-dimensional hy¬ 
drodynamic simulations have shown that the cooling phases that 
occur when these critical metallicities are reached can indeed trig¬ 
ger fragmentation in their associated density regimes (Bromm et al. 
(2001); Smith & Sigurdsson (2007); Smith et al. (2009); Meece 
et al. (2014) for the atomic and Clark et al. (2008); Dopcke et al. 
(2011, 2013) for the dust critical metallicity). However, it has also 
been found that the ability of gas to actually fragment depends 
sensitively on the initial conditions, no matter what the metallic¬ 
ity (Jappsen et al. 2007, 2009a,b; Meece et al. 2014). In particular, 
without a small amount of turbulence to create seed perturbations, 
monolithic collapse cannot be avoided (Meece et al. 2014). Thus, it 
seems that further progress requires understanding the true physical 
conditions in which the first metal-enriched stars formed. 

This paper is the first in a series that will investigate the range 
of scenarios in which star formation producing the near-universal 
IMF first occurred. Noting the ability of Pop III stars to completely 
evacuate their host mini-halos through ionizing radiation and su¬ 
pernovae (Whalen et al. 2004, 2008; Kitayama & Yoshida 2005), 
a number of recent works have sought to characterize the required 
conditions for metal-enriched star formation in terms of (i) incor¬ 
poration of metals via hierarchical structure formation into more 
massive halos (Greif et al. 2007; Wise & Abel 2008; Greif et al. 
2010; Wise et al. 2012b; Safranek-Shrader et al. 2014a,b, 2015) 
or (ii) fallback onto the original halo after after a recovery period 
(Ritter et al. 2012; Jeon et al. 2014; Ritter et al. 2014). These are 
both viable scenarios, but because they typically involve stars form¬ 
ing with metallicities of Z> 10“^ Z©, they are unlikely to be the 
formation mechanisms of the most metal-poor stars observed with 
[Fe/H]^ ~ — 5.5 (Norris et al. 2013). Here, we report on a third sce¬ 
nario for low-metallicity star formation, whereby mini-halos that 
have yet to form stars are externally enriched by a Pop III super¬ 
nova from a nearby halo. We note that while this channel has not 
before been shown conclusively to give rise Pop Il-like stars, it was 
previously suggested as a possibility by Greif et al. (2010), who 
found a halo similar to what we describe here but did not have the 
ability to follow the collapse and fragmentation. This is also similar 
to the “gravitational enrichment” scenario proposed by Maio et al. 
(2011). The simulations presented here are the first of their kind to 
follow the metal creation, non-uniform enrichment of the surround¬ 
ing medium, and eventual collapse to high density within a single. 


^ Unless stated otherwise, all densities and number densities are assumed 
to be in the proper frame. 

^ [Fe/H] refers to the log 10 of the iron abundance relative to the sun. 


coherent run. As such, they represent a significant leap forward in 
understanding the origins of the first low-mass stars. 

The paper is organized as follows. In Section 2, we describe 
the simulation code and the setup of the simulations. We present the 
results in Section 3. We first describe the general evolution of the 
simulation in Section 3.1. We then investigate the process by which 
the Pop II star-forming halo is enriched through turbulent mixing 
with the supernova blast-wave in Section 3.2. Following this, we 
characterize the physical conditions within the star-forming cloud, 
including fragmentation, the role of dust, and the velocity structure 
of the gas in Section 3.3. We conclude the results in Section 3.4 by 
estimating how often stars will form via this external enrichment 
mechanism. In Section 4, we discuss the implications of these re¬ 
sults and some limitations of our work. Finally, we summarize the 
results and discuss the ways in which this work may be extended in 
Section 5. 

2 METHODS 

2.1 Simulation Code 

2.1.1 The Enzo Framework 

We use the Enzo^ simulation code (Bryan et al. 2014) for all sim¬ 
ulations presented here. Enzo is an open-source, adaptive mesh 
refinement - 1 - N-body cosmological simulation code that has been 
heavily used for simulating cosmological structure formation over 
a large range of scales. Enzo has been employed by numerous 
works to study high-redshift structure formation, including Pop III 
star formation (Abel et al. 2002; O’Shea et al. 2005; O’Shea & 
Norman 2007, 2008; Turk et al. 2009, 2010, 2011a, 2012), low- 
metallicity star formation (Smith & Sigurdsson 2007; Smith et al. 
2009; Meece et al. 2014), and the first galaxies (Wise et al. 2012a,b; 
Xu et al. 2013; Wise et al. 2014; Chen et al. 2014). All functional¬ 
ity of the Enzo code is thoroughly detailed in Bryan et al. (2014), 
but we briefly describe the most relevant machinery here. Enzo 
uses a block-structured adaptive mesh refinement (AMR) frame¬ 
work (Berger & Colella 1989) to solve the equations of ideal hy¬ 
drodynamics in an Eulerian frame, with multiple hydrodynamic 
solvers implemented. Eor these simulations, we use the Piecewise 
Parabolic Method of Woodward & Colella (1984). This is coupled 
to an N-body adaptive particle-mesh gravity solver (Efstathiou et al. 
1985; Hockney & Eastwood 1988). When ionizing radiation fields 
from stellar sources are present (described below in Sections 2.1.2 
and 2.1.3), we also use the Moray adaptive ray-tracing radiation 
transport method of Wise & Abel (2011). The AMR framework 
dynamically creates (refines) and destroys (de-refines) grid patches 
at varying levels of resolution based on multiple criteria. This en¬ 
ables the study of astrophysical phenomena with processes occur¬ 
ring over a large range of spatial and temporal scales. We list the 
refinement criteria employed in this work below in Section 2.2. 

2.1.2 Chemistry and Cooling 

We use a chemistry model that follows the non-equilibrium evolu¬ 
tion of 12 primordial species (H, H^, He, He^, He^^, e“, H 2 , HJ, 
H“, D, D^, and HD) following the method originally described in 
Anninos et al. (1997) and Abel et al. (1997). This method has been 
updated to include H 2 formation via three-body reactions (Abel 

^ http://enzo-project.org/, changeset 7b5518a637 92 
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et al. 2002), H 2 formation heating, and collisionally-induced H 2 
emission (Turk et al. 2009). We include the effects of ionizing ra¬ 
diation from Pop III stars by coupling the chemistry network to the 
Moray radiation transport solver (Wise & Abel 2011). We also 
include optically thin, H 2 photo-dissociating Lyman-Werner (LW) 
radiation with an intensity decreasing as from the source to¬ 
gether with the H 2 “Sobolev-like” self-shielding model of Wolcott- 
Green et al. (2011). 

Just as primordial chemistry appears simple enough to allow 
for the hope that it can be computed entirely without limiting as¬ 
sumptions in a large, three-dimensional simulation, the introduc¬ 
tion of the first heavy elements by Pop III supernovae adds a level of 
complexity to the Universe that crushes that hope. We calculate the 
cooling from gas-phase metals following the method of Smith et al. 
(2008), which uses multidimensional tables pre-computed with the 
photoionization code, CLOUDY^ (Ferland et al. 2013). For sim¬ 
plicity, we assume the cooling from metals corresponds to a solar 
abundance pattern. When Pop III stars are present, we use a four¬ 
dimensional cooling table that spans the parameter space of density, 
metallicity, temperature, and electron fraction to account for the in¬ 
creased ionization from the radiation field. When metal-enriched 
gas begins to reach high densities, we switch to a table that covers 
only density, metallicity, and temperature, but was computed up to 
higher densities than the four-dimensional table. 

In addition to gas-phase chemistry and cooling, we also take 
into account the contribution of dust grains, which provide an addi¬ 
tional channel for H 2 formation and an efficient cooling mechanism 
at high densities. The implementation used here closely follows the 
work of Omukai (2000) and Omukai et al. (2005), with modifica¬ 
tions motivated by Dopcke et al. (2011), and is detailed in Meece 
et al. (2014). The primary assumption to note is that the gas-to-dust 
ratio scales with the metallicity and that the grain size distribution 
and abundances follow the model of Pollack et al. (1994), as in 
Omukai (2000). 


2.1.3 Population III Star Formation and Feedback 

We model the formation and feedback of Pop III stars with the 
method presented in Wise et al. (2012b). During the simulation, 
we insert a Pop III star particle into a grid cell where the following 
conditions are met: 

(i) The proper baryon number density exceeds 10^ cm“^. 

(ii) The gas fiow is convergent (V • Vgas < 0). 

(iii) The molecular hydrogen mass fraction, fH 2 = (ph 2 + 
PH 2 +)/Pb, exceeds 5 x 10“^, (where ph 2 , Ph 2 +’ Ph are the 
neutral molecular hydrogen, singly-ionized molecular hydrogen, 
and the total baryon densities, respectively). 

(iv) The metallicity is less than 10“® Z©. If the metallicity is 
greater than this, we will explicitly follow the collapse. 

These physical conditions are typical of metal-free molecu¬ 
lar clouds nearing the end of the ‘loitering’ phase observed in 
simulations of Pop III star formation (Abel et al. 2002; Bromm 
et al. 2002; Yoshida et al. 2006; O’Shea & Norman 2007; Turk 
et al. 2009). In the simulations here, we adopt a delta function 
at M = 40 Mq for the IMF of the Pop III star particles. When 
a star particle is created, the mass of the star is removed evenly 
from a spherical region containing twice the particle mass. Once 
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formed, the star particle isotropically emits ionizing and LW ra¬ 
diation following the stellar evolution models of Schaerer (2002) 
with no mass loss. Following this, the main-sequence lifetime of 
the 40 Mq Pop III star is 3.86 Myr. Furthermore, we treat the 
star as a radiation point source, whose spectrum is discretised into 
four energy groups with the following photon energies and lumi¬ 
nosities: (i) hydrogen ionizing photons with = 28 eV and 
= 2.47 X 10^^ s (ii) helium singly-ionizing photons with 
= 30 eV and — 1.32 x 10^^ s“^; (iii) helium doubly- 
ionizing photons with E' = 58 eV and = 8.80 x 10^® s“^; 
(iv) H 2 dissociating photons with — 2.90 x 10^^ s“^ (Schaerer 
2002). The ionizing radiation field is transported using the MORAY 
radiation transport solver (Wise & Abel 2011). As stated previ¬ 
ously, we model the H 2 photo-dissociating LW flux to be declining 
simply as r“^ with a self-shielding correction. 

When the Pop III star’s main-sequence lifetime ends, it ex¬ 
plodes as a standard Type II core-collapse supernova with an ex¬ 
plosion energy of 10^^ ergs. Following the calculations of Nomoto 
et al. (2006), the supernova has a metal yield of 11.19 M© and a 
total ejecta mass of 38.6 M©. We model the blast-wave by deposit¬ 
ing the explosion energy and mass into a sphere with a radius of 
10 proper pc that has been smoothed at its surface to improve nu¬ 
merical stability (Wise & Abel 2008). We choose this radius as it 
marks the end of the free-expansion phase of the supernova. At the 
time of instantiation, the blast-wave is extremely well-resolved and 
agrees with the Sedov-Taylor solution. 

2.2 Simulation Setup 

We initialize the simulation with a 500 kpc/h comoving box at 
z = 180 using the MUSIC initial conditions generator (Hahn & 
Abel 2011) with the WMAP 7 best-fit cosmological parameters, 
FLrri = 0.266, Fix = 0.732, = 0.0449, Ho = 71.0 km/s/Mpc, 

as — 0.801, and n^ = 0.963 (Komatsu et al. 2011), and us¬ 
ing a Eisenstein & Hu (1999) transfer function and second-order 
Lagrangian perturbation theory. We first run an exploratory dark 
matter-only simulation with 512^ particles to 2 : = 10. We then 
use the ROCKSTAR halo finder (Behroozi et al. 2013a, inter¬ 
faced through the YT analysis toolkit) and the CONSISTENT-TREES 
merger tree algorithm (Behroozi et al. 2013b) to locate a halo with 
total mass (at z = 10) of ~ 10^ M© and the highest number of 
unique progenitor halos with M ^ 10® M©. The second of these 
criteria was added to maximize the number of Pop III stars formed, 
and hence the amount of metals produced, by the time the halo has 
reached z = 10. The halo selected, referred to from here on as the 
“target halo”, has a total mass dX z — 10 of 1.7 x 10^ M© and 4 
unique progenitors with M > 10 ® Mq (and 7 with M ^ 5 x 10® 
Mq.) 

We generate the initial conditions for the primary simulation 
using the same configuration and random seed as the exploratory 
run, including baryons as well as dark matter, and with two addi¬ 
tional levels of telescoping refinement around the target halo. The 
high resolution region is a rectangular prism minimally containing 
all of the dark matter particles that end up within 3 virial radii of the 
target halo. It has an effective resolution of 2048®, corresponding 
to a comoving spatial resolution of 0.244 kpc/h, and a baryon (dark 
matter) mass resolution of 0.259 M© (1.274 M©). 

We evolve the simulation from z = 180 with the physics de¬ 
scribed in Section 2.1, allowing adaptive mesh refinement to occur 
only within the high resolution region. Grid cells are adaptively re¬ 
fined by factors of 2 in each dimension when any of the following 
conditions exist: 
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(i) The dark matter mass within a grid cell is greater than 4 times 
the initial mass (i.e., when more than 4 dark matter particles are in 
the same cell). 

(ii) The gas mass within a grid cell is greater than 4 times the 
mean baryon mass per cell on the root grid multiplied by a factor, 

where L is the refinement level. 

(iii) The local Jeans length is resolved by less than 64 cells. 

The negative exponential factor used with the baryon mass criterion 
ensures that refinement is mildly super-Lagrangian. Requiring that 
the local Jeans length be resolved by at least 64 cells (a factor of 16 
higher than the criterion put forth by Truelove et al. 1997) has been 
shown to be necessary to guarantee that fragmentation is purely 
physical in origin (Federrath et al. 2011; Meece et al. 2014). 

We place no strict limit on the maximum number of AMR lev¬ 
els in the simulation. However, in practice, the simulation typically 
reaches 15 levels of refinement before forming a Pop III star, corre¬ 
sponding to a comoving spatial resolution of 0.03 pc/h. Finally, we 
stop the simulation once gas with metallicity of at least 10“® Z© 
reaches a maximum number density of ~ 10^^ cm“^. This occurs 
at 27 levels of refinement, for a spatial resolution of approximately 
1.5 AU/h comoving. As discussed in Section 3.3, after the first Pop 
III supernova occurs, we start a second simulation from this point 
with all physics the same except without dust grains. 

We run the simulations on the Blue Waters supercomputer at 
the National Center for Supercomputing Applications on 64 nodes 
of the machine, using 8 cores/node. The primary run took roughly 
42 days, with a total computational cost of approximately 64,000 
node-hours. For all of the analysis and images presented here, we 
use the YT^ simulation analysis toolkit (Turk et al. 2011b). 

3 RESULTS 

3.1 Qualitative Evolution 

The two simulations proceed identically until the collapse of metal- 
enriched gas first occurs at z ~ 16.6, at which time their behaviors 
diverge due to the presence/absense of dust grains. The qualita¬ 
tive evolution prior to this is as follows. From here on, we refer 
to the halo in which metal-enriched collapse occurs as the “action 
halo.” Within the high resolution region, two Pop III stars form in 
separate halos, shine, and explode. The evolution of the halos host¬ 
ing Pop III stars agrees well with earlier studies of this scenario 
(Whalen et al. 2004, 2008; Kitayama & Yoshida 2005). As the ha¬ 
los are only a few hundred thousand M©, the Pop III star is able 
to completely photo-evaporate them. The central baryon density is 
reduced to roughly 0.1 cm“^ before the core-collapse supernova 
occurs and totally evacuates the gas from the halo. As in Whalen 
et al. (2008), the supernova shockwave has a density of just under 
1 cm“^ and catches up with the expanding shell of the HII region 
within a few million years of the explosion. This collision enhances 
the density by a factor of roughly 10 to 30, after which the shell pro¬ 
ceeds to break up into small knots. These knots eventually dissolve, 
but they carry metals out into the intergalactic medium (IGM), cre¬ 
ating finger-like structures in the metal field, which are clearly vis¬ 
ible in the upper-right panel of Figure 13. The first of the two Pop 
III stars, forming at z ~ 23.7, has no influence on the action halo. 
The second Pop III star forms at z ~ 18.2 (~ 66.6 Myr later) and is 
solely responsible for enriching the action halo. By the time of gas 

^ http://yt-project.org 


collapse in the action halo, the supernova blast-waves from each of 
the two Pop Ill-hosting halos (separated by ~ 3.4 kpc) are still ~ 
1 kpc from a collision. In another version of this simulation carried 
beyond this point, we find that the blast-waves meet at z ~ 14 (~ 
65 Myr in the future). 

Figure 1 displays the sequence of events immediately leading 
up to metal-enriched collapse. The action halo is roughly 200 pc 
away from the second Pop III star when it forms. The Pop III star 
photo-ionizes its host halo and causes a slight compression of the 
gas in the filament hosting the action halo. This is visible in the third 
row of Figure 1 . Prior to this, the action halo has a relatively low 
central baryon fraction, about a few percent of the cosmic baryon 
to dark matter ratio. Over the main sequence lifetime of the star, 
the central density of the action halo increases by a factor of a few, 
raising the baryon fraction to over 10% of Qb/^dm- The gas in the 
action halo remains almost entirely neutral, with an HII fraction 
of about 10“^ in the central 10 pc, increasing to 10“^ at 50 pc. 
Approximately 6 Myr after the Pop III supernova occurs, the blast- 
wave encounters the action halo moving in the opposite direction. 
At this point, the action halo has a virial mass of ~ 3 x 10^ M©. 
Turbulence produced by the virialization of the action halo (Wise 
& Abel 2007; Greif et al. 2008) and the passing of the blast-wave 
allows metals from the supernova to mix into the core of the action 
halo. It also happens that the blast-wave encounters the action halo 
just prior to colliding with the HII region shell and thus meets the 
halo with an enhanced density. This process is illustrated in Figure 
2. Finally, runaway collapse occurs in the center of the action halo 
where the metallicity has reached ~ 2 x 10-® Zq. While this is 
happening, the action halo and the gas-depleted Pop III halo con¬ 
tinue on a near collision course. By the time collapse occurs, the 
two halos have passed close enough that the halo finder considers 
them to be a single halo, although distinct cores are still visible. 

3.2 Enrichment of the Action Halo 

The Pop III supernova blast-wave carries with it 11.19 M© of met¬ 
als and impacts the action halo about 6 Myr after the explosion. 
Figure 3 shows the metallicity as a function of radius from the 
metal-enriched gas cloud at the time of final collapse. At the scale 
of 10-100 pc, gas exists from zero metallicity up to roughly 10 Z© 
as spherical shells of this size are sampling the blast wave itself as 
well as regions that have yet to be encountered. Within a few pc of 
the dense cloud, no pristine gas exists and the metallicity is as high 
as 10“^ Z©. Inside of 0.01 pc, the metallicity is nearly constant 
at -2x10“^ Z©. To understand this, we compare the timescales 
for collapse and mixing in Figure 4. The collapse timescale is the 
free-fall time, given by 

^ \l32Gp’ 

where p is the gas density. We estimate the mixing timescale as the 
vortical time (e.g., Ritter et al. 2014), 

“ W^\’ 

where the factor of 27r accounts for a full rotation of an eddy. The 
mixing timescale is shorter than the free-fall time outside of a few 
pc and, thus, mixing occurs here. However, at r ~ 1 pc, the free-fall 
time is slightly shorter than the vortical time, allowing the gas to 
momentarily collapse faster than the rate at which additional metals 
can be brought in from the surrounding medium. We expect the rate 
of collapse to be regulated by the free-fall time in this regime since 
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Figure 1. Slices illustrating the large-scale evolution of the simulation prior to the collapse of metal-enriched gas, oriented such that it contains the vector 
connecting the Pop III star to the rough center of the action halo. Letters indicate the following: at 2 ; = 19.287 a) the Pop III star-forming halo, b) the halo in 
which metal-enriched gas will collapse (i.e., the action halo), c) the Pop III star begins shining ~ 17 Myr later, d) the Pop III star explodes in a core-collapse 
supernova ~ 4 Myr later, and e) the supernova blastwave collides with the action halo. The black circles in the top-left panel indicate the virial radii of the 
Pop III (79 pc) and action (82 pc) halos. The white circle in the bottom-left panel shows the location of the evacuated Pop III star halo. 
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Figure 2. Slices showing the turbulent mixing of the metal-enriched blastwave with the action halo. The top panels show the state of the simulation ~ 18.75 
Myr after the bottom panel of Figure 1. Here, as in Figure 1, the scale bars show distances in the proper frame and the slicing plane contains both the Pop III 
star and action halo center. 
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Figure 3. Phase plot showing the mass of gas in bins of metallicity and ra¬ 
dius from the central core for the simulation with dust. All gas with metal¬ 
licity Z < 10“^ Zq is combined and shown in the small strip at the bottom. 
This gas is entirely at log(r/pc) > 0.5; all gas within that radius has been 
enriched to a higher level. 
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Figure 5. The red line shows the maximum metal mass that would have 
passed through a circle of radius, r, in the absence of mixing after the blast- 
wave has reached the action halo. The green line shows the enclosed mass 
of metals within the cloud and blue shows the enclosed metal mass corre¬ 
sponding to a constant metallicity of 2 x 10“^ Zq. The flattening of the 
red line indicates the total metal yield of the supernova (11.19 Mq). All of 
the metals are within 500 pc of the center, as shown by the flattening of the 
green line. 


the cooling time is equivalent to or shorter than the free-fall time 
here. 

Further evidence of the mixing process can be seen in Figure 
5, where we compare the enclosed mass in metals to the maximum 
flux of metals carried by the blast-wave. Under the assumptions that 
the blast-wave is spherical and that all of the metals are carried in 
the front, the mass in metals that would have passed through a hole 
of radius, r, is given by 
^2 

Mflux = Msn 7 ^, (3) 

where the metal yield of the supernova, Msn, is 11.19 Mq, and the 
proper distance between the halo and the origin of the blast-wave. 


d, is ~70 pc. Where the enclosed metal mass is greater than the 
amount carried by the blast-wave, turbulence has acted to enhance 
the metal mass over what the blast-wave could have maximally de¬ 
posited without mixing. We also show in Figure 5 the mass in met¬ 
als corresponding to the constant metallicity of the inner cloud. In 
the region where this curve diverges from the true enclosed metal 
mass, mixing is ongoing and continues to increase the metallicity 
of the gas beyond the point at which the metallicity is frozen out by 
rapid collapse. We note here for the sake of novelty that the total 
mass in metals within 1 pc of the center of the dense cloud is only 
about half the mass of Jupiter. 

Cen & Riquelme (2008) study the mixing of metals into pris¬ 
tine halos impacted by supernova blast-waves using idealized sim¬ 
ulations. In the conflguration most similar to the scenario here, 
(Mhaio ~ 3 X 10^ Mq, V ~ 11 km/s,) they And that very little 
mixing occurs with a majority of the dense gas in the halo remain¬ 
ing at a metallicity less than 3% of the blast-wave. However, there 
are some key differences between the conditions we observe in the 
action halo and the simulations of Cen & Riquelme (2008). The 
minimum halo mass simulated in Cen & Riquelme (2008) was 10® 
Mq, roughly 3 times more massive than the action halo. More im¬ 
portantly, the halos in Cen & Riquelme (2008) were initialized with 
the gas in hydrostatic equilibrium with no internal motion. As Fig¬ 
ures 1 and 2 show, the action halo is in a highly perturbed state 
immediately prior to impact, having its turbulence first stirred dur¬ 
ing virialization and then being eroded by the Pop III stellar radi¬ 
ation. This is supported by the results of Richardson et al. (2013), 
who perform a similar study to that of Cen & Riquelme (2008), 
but using realistic mini-halos with M ~ 10® Mq extracted from 
a cosmological simulation. Whereas Cen & Riquelme (2008) And 
that most of the material in idealized halos of this size remains pris¬ 
tine, Richardson et al. (2013) observe that most of the halo gas is 
enriched to ~ 3% of the metallicity of the blast-wave. Finally, we 
note that while enough metals are mixed into the collapsing gas to 
constitute a transition to low-mass star formation in the cloud, the 
mixing should not be considered very efficient. At the 1 pc scale 
where free-fall collapse freezes in the metallicity of the dense gas. 
Figure 3 shows that the range in metallicity is nearly 3 orders of 
magnitude, with metallicities present up to 10 times higher than 
that of the dense cloud. 


3.3 Physical Conditions of The First Pop II Stars 

The first Pop II star forming cloud within the simulation has a 
nearly constant metallicity of Z ~ 2x10“® Zq. In Figure 6, we 
show density projections progressively zooming in to the central 
density peak of the final output of the simulation. Not surprisingly, 
we see no evidence of fragmentation at densities much lower than 
the maximum reached in the simulation, as gas-phase metal cool¬ 
ing has almost no effect at such low metallicities (Omukai 2000; 
Omukai et al. 2005; Smith et al. 2009; Meece et al. 2014). How¬ 
ever, in the central 100 AU, we observe significant fragmentation at 
densities, p > 10“^^ g/cm®, where the cooling from dust grains be¬ 
comes efficient (Omukai 2000; Omukai et al. 2005; Schneider et al. 
2006; Clark et al. 2008; Dopcke et al. 2011; Meece et al. 2014). To 
confirm that the fragmentation that occurs is due to the effects of 
dust, we run a second simulation without dust starting from the 
moment following the first Pop III supernova. In the lower-right 
panel of Figure 6, we show the corresponding central, dense region, 
which has undergone significantly less fragmentation. 

To quantify this notion, we run a clump-finding algorithm to 
identify bound clumps within 50 pc of the peak density in both sim- 
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Figure 4. Top: radial profiles from the final output of the simulation with dust of the mass-weighted mean free-fall time (red), cooling time (green), and vortical 
time (blue). Bottom: ratio of cooling and vortical times to the free-fall time. 


ulations. Clump identification (Smith et al. 2009; Turk et al. 201 lb) 
proceeds by using a contouring algorithm to calculate the largest 
topologically disconnected structures that satisfy the following cri¬ 
terion: 

KE + TE - y] Ai idyn.i < PE, (4) 

i 

where KE, TE, and PE are the total kinetic, thermal, and potential 
energies of the clump, and Ai and tdyn,i are the cooling rate and 
dynamical time for each member grid cell (Meece et al. 2014). As 
in Meece et al. (2014), we include the cooling term in Equation 4 
to account for marginally unbound clumps that are rapidly cooling 
and, hence, expected to become gravitationally bound in the near 
future. In Eigure 7, we show the number of clumps as a function 
of the minimum gas density within the clump. The minimum den¬ 
sity within the clump effectively serves as the density at which the 
fragmentation event forming that clump occurred. In the simula¬ 
tion with dust, we identify 126 individual clumps, compared with 
38 in the simulation without dust. Of the 38 clumps in the simula¬ 
tion without dust, 23 are at densities less than 10“^^ g/cm^, similar 
to the simulation with dust, which has 21 in the same range. The 
masses of these are also similar, ranging from about 0.1 M© up to 
about 1000 M©. The high density clumps (pmin > 10“^^ g/cm^) 
in the simulation with dust show a roughly log-normal distribution 
in masses, centered at ~ 10“^ M© and ranging from about 10“® 
M© to about 10“^ M©. The total gas mass within the central 50 AU 
is about 0.8 M©. While is is likely that a number of these may later 
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Figure 7. The number of bound clumps at a given threshold density (i.e., the 
minimum density of the clump) within the central 50 pc of each simulation. 



merge, it is clear that the presence of dust significantly increases 
the degree of fragmentation during this phase of the collapse. 

The increased fragmenation observed here can be conclusively 
attributed to the presence of dust by comparing the thermal state 
between between the simulations with and without dust, as we 
do in Figure 8. The sharp decrease in temperature seen at high 
density, n > 10^^ cm“^, occurs when the dust and gas become 
thermally coupled and continuum emission from the dust is able 
to efficiently cool the gas. With dust present, gas in the density 
range, 10^ cm“^ < n < 10^^ cm“^, is also cooler by a factor of 
roughly 2 owing to the increased H 2 fraction from the availability 
of the dust-grain formation channel, as shown in Figure 9. These 
results are all in reasonable agreement with analagous one-zone 
models (Omukai 2000; Omukai et al. 2005; Schneider et al. 2006) 
and three-dimensional simulations (Clark et al. 2008; Dopcke et al. 
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Figure 6. Mass-weighted density projections showing the collapse of the metal-enriched molecular cloud from the final output of the simulations. The bottom 
panels show the central, dense core of the simulation with (left) and without (right) dust. Distances shown are in the proper frame. 


2011; Meece et al. 2014). The simulation without dust also shows 
a decrease in temperature at high density, albeit much less pro¬ 
nounced. As Figure 9 shows, this is due to the rapid increase in H 2 
fraction caused by 3-body reactions becoming efficient. We note 
that this appears to happen somewhat later than in previous works 
(e.g., Omukai et al. 2005). However, it was shown by Meece et al. 
(2014) that the precise timing of collapse can significantly affect 
chemical evolution within the gas. 

Finally, we briefiy examine the velocity structure of the col¬ 
lapsing gas cloud in Figure 10. Here, as in Figures 4 and 5, we 
calculate the center of the cloud by creating a sphere of 100 pc cen¬ 
tered on the halo, and iteratively recalculating the center of mass 
after decreasing the radius of the sphere by 10% until reaching a 
radius of 1 AU proper. For the velocity profiles, we subtract off the 
bulk velocity of the central 1 AU. We define the rotational velocity 
as the component that is orthogonal to both the radial and angular 
momentum vectors, where the angular momentum vector is calcu¬ 
lated from the central 1 AU. The cloud is dominated by turbulent, 
or at least unordered, motion as the variance of the velocity mag¬ 
nitude (minus the bulk velocity of the center) is higher than both 
the radial and rotational components at all radii. At large radii, the 
velocity variance is between 1.5 and 2 times the sound speed and 
closely follows the sound speed at lower radii, only exceeding it at 
r ~ 10 AU where the cooling leads to rapid cooling. Though less 
than the variance at all radii, there is still a significant infall com¬ 
ponent that is generally trans-sonic. We find that in the simulation 


without dust where the gas temperature slightly higher, the infall 
is similarly raised. The gas shows almost no coherent rotation with 
an average rotational velocity that is near zero at virtually all radii. 
This is not simply an effect of where we have chosen to calculate 
the angular momentum vector as we find that calculating it at r = 
10 and 100 AU gives qualitatively similar results. 


3.4 The Bigger Picture for Metal-Enriched Star Formation 

The gravitational binding energy of mini-halos is well below the 
energy of even a standard core-collapse supernova. As such, it 
will be some time before metal-enriched stars are able to form in- 
situ from self-enrichment of the progenitor halos that form Pop III 
stars. However, just two Pop III stars have enriched a considerable 
amount of space and material by z ~ 16.6, when this simulation 
stops. In Figures 11 and 12, we show the volume and mass en¬ 
riched to metallicities high enough to transition to a metal-enriched 
star formation mode. Roughly 1% of the mass and just under 1% 
of the volume within the high-resolution region of the simulation 
are enriched to a metallicity of 10“^ Z© by the final output. Only 
a small fraction of the metals created in the two supernovae will 
be captured by nearby halos and be incorporated into stars by the 
time the target halo has assembled at z = 10 with a total mass of 
~2 X 10^ M©, about 250 Myr from the current time in the simula¬ 
tion. If metal-enriched star formation occurs on free-fall timescales, 
this 250 Myr time period corresponds to a density of roughly 10“^^ 
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Figure 10. Radial profiles of the mass-weighted, mean velocity magnitude (red), rotational velocity (green), radial velocity (blue), and sound speed (purple). 
The variance of the total velocity is shown in grey and the circular velocity corresponding to the enclosed mass is shown in black. 


glow?. Following this assumption, we can estimate the number of 
additional halos externally enriched by these two Pop III super¬ 
novae that will be able to form Pop II stars in this time period by 
searching for all the halos with metallicities of at least Z© in 
gas with 10 “^^ g/cm^. 

In Figure 13, we show the results of this attempt to predict 
the future. In this exercise, we consider halos with masses as low 
as 10^ M 0 . We identify 21 halos in which Pop II stars could form 
before the target halo has assembled. In cases where the metallicity 
is less than ~ 10 “^ Z©, gas-phase metal cooling will be insuffi¬ 
cient at low densities, requiring a critical H 2 fraction in a scenario 
analogous to the filtering mass for Pop III star formation (Tegmark 
et al. 1997; Yoshida et al. 2003). To account for this, we require 
halos with Z < 10~^ Z© to have a minimum mass of 10^ M©. In 
other words, halos with M < 10^ M© must have metallicities of at 
least 10“^ Z 0 such that they are able to cool via metal lines. This 
constraint lowers the number of candidate halos to 16. Of these 16 
halos, 7 have masses of 10^ M 0 < M < 10^ M 0 , 6 have masses 
of 10^ M 0 < M < 10^ M 0 , and 3 have masses of 10^ M 0 < 
M < 10® M 0 . Despite having a 66 Myr head start, only 1 of the 
16 candidate halos exists in the region enriched by the first of the 
two supernovae. In fact, this halo is the original host of the Pop III 
star. Out of random chance, the halo hosting the first Pop III star 
has significantly fewer neighbors than that of the second Pop III 
star. It has been approximately 96 Myr since the explosion of the 
first supernova. That this halo has not yet recovered is consistent 


with Jeon et al. (2014), who find a recovery time of 140 Myr for a 
3x10® M 0 halo experiencing a 10®^ erg explosion. The lower pan¬ 
els of Figure 13 show that the halo hosting the second Pop III star 
(as well as the action halo) formed in a significantly denser region 
with many more neighboring halos. Out of the 22.39 M 0 of metals 
created, about 1.8 M 0 have been captured by neighboring halos. It 
is important to note that as metal-enriched stars form within region 
b of Figure 13, their feedback may infiuence other nearby candidate 
halos, either inhibiting star formation with radiation or contributing 
additional metals via supernovae. Thus, the resulting star formation 
in this region may be slightly different from the picture presented 
here, although it is difficult to know for sure without running the 
simulation into the future including Pop II star particles. 

What is most striking is the range in metallicities resulting 
from a single enrichment event. A total of 7 halos have metallici¬ 
ties in their dense gas (p ^ 10 “^® g/cm®) of at least 10 “® Z 0 , with 
a single halo as high as 0.1 Z 0 . The halo with the highest metal¬ 
licity is, in fact, the combined action/Pop III halo, represented in 
Figure 13 as the large, red circle in the center of region b. The core 
of the Pop III halo is just barely visible to the lower-right of the 
main core of this structure. Just after the collision with the blast- 
wave but before encountering the Pop III halo, the action halo ex¬ 
periences a major merger. The final object, including the not yet 
totally merged Pop III halo, has a mass of 9.8 x 10® M 0 . Interest¬ 
ingly, it seems that it will likely form both the lowest and highest 
metallicity Pop II stars from the metals created thus far. More im- 
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Figure 13. Top panels show projections of overdensity (left) and metallicity (right) in a region minimally containing all of the metals produced in the simulation. 
Bottom panels show projections of the regions surrounding each of the two supernova events. Region a contains the first Pop III star to form and region b 
contains the second as well as the action halo. The locations of halos with metallicities of at least 10“® Zq in dense gas (p ^ 10“^^ g/cm^.) The size and 
color of the circles indicate the virial radii and peak dense gas metallicities, respectively. The dashed lines indicate halos with masses less than 10^ Mq and 
metallicities below 10“"^ Zq which may be unable to cool. 


portantly, the picture presented here shows that the first generation 
of metal-enriched stars was highly heterogeneous in metallicity, en¬ 
compassing significantly more than just the most metal-poor stars. 


4 DISCUSSION 
4.1 Implications 

The one-zone collapse models with their extremely sophisticated 
chemistry networks (Omukai 2000; Omukai et al. 2005, and re¬ 
lated works) have paved the way toward understanding the required 
conditions for fragmentation from the perspective of the ability of 
low metallicity gas to cool. The simplicity of the collapse model 
has enabled the exploration of a large parameter space of non-solar 
abundance patterns motivated by Pop III supernova yield predic¬ 
tions (Schneider et al. 2006) and varying levels of dust (Schneider 
et al. 2012). Encouragingly, the idealized three-dimensional simu¬ 
lations (Bromm et al. 2001; Smith & Sigurdsson 2007; Clark et al. 
2008; Smith et al. 2009; Dopcke et al. 2011, 2013; Meece et al. 
2014) have confirmed that the thermal evolution predicted by these 
models can, indeed, induce fragmentation provided there is a small 


degree of turbulence to seed density perturbations. This is particu¬ 
larly interesting given that the recently discovered most metal-poor 
star, SMSS 1031300.36-670839.3 (Keller et al. 2014), has abun¬ 
dance patterns compatible with the yield of a single Pop III super¬ 
nova (Keller et al. 2014; Maras si et al. 2014). Additionally, Maras si 
et al. (2014) have shown that, under a range of dust abundances, 
a collapsing gas-cloud with the chemical composition of SMSS 
J031300.36-670839.3 would have undergone low-mass fragmenta¬ 
tion. What has been missing from this picture is a means to implant 
such low metallicities into a star forming environment. We have 
provided such a mechanism in this work. 

The frequency with which this mechanism operates is uncer¬ 
tain, although it is also encouraging that Greif et al. (2010) appear 
to have identified a similar situation, albeit with a higher metallicity 
(Z-2xl0-^ Zq). We have, so far, only followed this one occur¬ 
rence, although we identify a number of other externally enriched 
halos that could potentially form stars with metals from a single 
supernova. However, it is clear that the number of halos externally 
enriched by a given Pop III supernova will vary greatly. For exam¬ 
ple, the region around the first Pop III star in our simulation has 
no viable externally enriched halo candidates and will see second 
generation star formation only after the Pop III halo itself recovers. 
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Figure 8. Phase plots of the gas mass in bins of number density and tem¬ 
perature for the simulation with (top) and without (bottom) dust. In the top 
panel, the red shows the analogous distribution for the dust temperature. 
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Figure 9. Prohles of the mass-weighted average H 2 fraction as a function 
of number density for the two simulations. 


By contrast, the region around the second Pop III star has a num¬ 
ber of such candidates. To first order, the likelihood of the external 
enrichment mechanism depends on the average separation of Pop 
III star-forming halos from mini-halos that have yet to undergo star 
formation. We provide a crude estimate of this in Figure 14. Us¬ 
ing the exploratory simulation, we identify all halos with masses 
greater than 10^ M© and calculate the distance to the nearest halo 
with a mass in the range of 10 ^ M© ^ M ^ 5 X 10® M©. The 


z 



t [Myr] 


Figure 11. The volume within the simulation enriched to at least a given 
metallicity as a function of redshift. The different colors show metallicities 
ranging from Z = 10“® Zq (blue) to 10 Zq (green) increasing by factors of 
ten. The red line shows 1% of the total volume of the high resolution region. 
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Figure 12. The mass of gas enriched to at least a given metallicity as a 
function of redshift. The different colors show metallicities ranging from Z 
= 10“® Z 0 (blue) to 10 Zq (green) increasing by factors of ten. The red 
line shows 1 % of the total baryon mass within the high resolution region. 
Note, the refinement region is in a large-scale overdensity and has roughly 
2.8 times more mass than an average region of that size. 


former group is assumed to have a high probability of forming a 
Pop III star, where the latter may be just about to do so. We allow 
the mass ranges to overlap noting that the action halo is actually 
more massive than the Pop III halo. During the period from z = 20 
to z = 10, the average separation of these pairs is roughly 500 pc 
and evolves very little. The high resolution region of the simulation 
shows a similar evolution but with a slightly lower value of about 
300-400 pc owing to the fact that it is located within a large-scale 
overdensity. The distance between the Pop III halo and the action 
halo is slightly more than 1 standard deviation lower than the aver¬ 
age pair separation within the full simulation volume. This suggests 
that the external enrichment mechanism is a non-negligible compo¬ 
nent of low-metallicity star formation. Jeon et al. (2014) find that 
the recovery time for mini-halos with M < 5 x 10^ Mq (similar 
to the Pop III star-forming halos in our simulations) is of the order 
of 100 Myr. In the simulations presented here, the time between 
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Figure 14. The blue line and shaded region shows the average distance be¬ 
tween halos with mass, M ^ 10^ Mq , and the nearest halo with mass, 10^ 
Mq ^ M ^ 5 X 10^ Mq as a function of redshift for the entire simulation 
volume. The red line and shaded region shows the measurement for just the 
high resolution region. The green line shows the distance between the ac¬ 
tion halo and the Pop III star-forming halo. The three black stars indicate 
the moments when the Pop III star forms, explodes, and when the supernova 
blast-wave collides with the action halo. 

the explosion of the second supernova and the collapse of metal- 
enriched gas is roughly 25 Myr. The shorter timescale is not sur¬ 
prising, as the mixing of metals into a nearby halo will occur early 
on when the blast-wave is still moving outward quickly. Thus, if 
the conditions are right for the external enrichment mechanism to 
occur, it should typically produce metal-enriched stars before they 
can form from fallback of metals onto the Pop III halo. 

Finally, the candidate halos surrounding the second Pop III 
star have a wide range of metallicities, from just over 10“^ Zq to 
nearly 10“^ Zq. This range of metallicities is, itself, a striking re¬ 
sult as it implies that there may be stars that are only moderately 
metal-poor by observational standards (with, e.g., thousands exist¬ 
ing in the SDSS SEGUE sample) yet trace a single progenitor. This 
could present a new opportunity to gain access to the Pop III IMP, 
although one that will certainly not be without its own challenges. 

4.2 Caveats and Limitations 

We make strong assumptions about the properties of the metals and 
dust grains in these simulations. We assume that gas-phase met¬ 
als have a solar abundance pattern and that the dust grains are of 
a composition, size distribution, and fraction of total metal mass 
identical to the interstellar medium of the Milky Way. While these 
assumptions are not explicitly correct, the thermal evolution we ob¬ 
serve is still in reasonable agreement with studies that use realistic 
abundance patterns and dust properties (e.g. Schneider et al. 2006, 
2012; Marassi et al. 2014). We also assume that the dust grain pop¬ 
ulation is static after the inital grain formation, although it has been 
shown recently that grains can continue to grow during collapse 
and that their effects can be important even when the initial deple¬ 
tion fraction is very small (Nozawa et al. 2012; Chiaki et al. 2013, 
2014). These effects, along with more realistic abundance patterns, 
will be important to consider in the future when these simulations 
can be carried to the point where an IMP can be reasonably ex¬ 
trapolated from their final state. However, given the current lack 
of understanding of the physical origin of the stellar IMP, such an 


exercise is not warranted at this time. In relation to the chemistry 
model employed here, the primary goal of this work is to show 
that the high density cooling phase associated with dust can induce 
fragmentation in realistic conditions. 

We have chosen the formation of a 40 Mq Pop III star fol¬ 
lowed by a 10^^ erg core-collapse supernova as the sole outcome 
of the collapse of metal-free gas. Since this work deals with enrich¬ 
ment of star forming gas by a single supernova, this assumption is 
as valid as choosing from a realistic IMP (e.g. Hirano et al. 2015) 
provided that the mass chosen is plausible within that IMP, which 
in this case is true. However, given the range of alternative Pop III 
explosion scenarios, such as energetic hypernova (Nomoto et al. 
2006) and pair-instability supernova (Heger & Woosley 2002), de¬ 
termining the robustness of the external enrichment mechanism to 
different explosion energies and yields has significant merit and 
is a topic we will explore in a future study. Por example, Cooke 
& Madau (2014) claim that the lack of observational evidence for 
Pop III pair-instability supernova does not necessarily rule out their 
existence entirely. Instead, they argue that the transition to Pop 
II star formation may simply be dominated by low energy super¬ 
novae whose ejecta fall back quickly while pair-instability events 
are masked by the fact that they evacuate their host halos and have 
very long recovery times. However, if pair-instability supernovae 
can just as easily give rise to the external enrichment mechanism 
presented here, it could present a problem for this explanation. 

Given the nature of our code, which is both grid-based and 
uses adaptive mesh refinement, it is possible that the mixing of 
metals is artificially enhanced in our second-generation halo. This 
is a known property of grid-based codes (e.g., Wadsley et al. 2008; 
Joggerst et al. 2014), and must be kept in mind when interpreting 
results such as those presented in Pigure 3. In this particular sim¬ 
ulation, however, there is reason to believe that numerical mixing 
is not the origin of the uniform metallicity seen in the center of 
our second-generation halo. This is due to the nature of the fiuid 
fiow in the halo, which is highly turbulent. Turbulence results in 
the efficient mixing of metals over the length scales of the tur¬ 
bulent cascade on the scale of a few local eddy turnover times, 
which are comparable to the local dynamical time in supersonic 
turbulence (Balsara & Kim 2005; Pan & Scannapieco 2010). As 
can be seen in Pigure 4, the dynamical, cooling, and vortical (i.e., 
turbulence) scales are quite comparable over the inner parsec of 
the second-generation halo, suggesting that the gas is well-mixed. 
This thorough mixing may also be observed in nature - Karlsson 
et al. (2012) shows the tentative existence of a dissolved star cluster 
with [Pe/H] ~ —2.7 in the Sextans dSph galaxy. This star cluster 
was found using chemical tagging, and found a handful of very 
old stars that are highly clustered in a multi-dimensional chemical 
abundance space. This clumping indicates chemical homogeneity 
in the gas cloud out of which these stars formed, which is estimated 
to be substantially more massive than the one under consideration 
in our calculation, and supports the assertion that gas should be 
well-mixed due to turbulent processes. 


5 SUMMARY AND CONCLUSIONS 

We present the first results from a study designed to characterize 
the physical conditions associated with the transition from Pop III 
to Pop II star formation. We have simulated the formation of the 
first Pop II stars in the Universe in a single, coherent simulation be¬ 
ginning from cosmological initial conditions. The simulations have 
a dark matter particle mass of just under 1.3 Mq, allowing us to 
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well-resolve the small (few x 10^ M©) mini-halos that host the first 
stars. We insert 40 M© Pop III star particles in metal-free molec¬ 
ular clouds that have evolved past the loitering phase and simulate 
their main-sequence lives with radiation hydrodynamics. Each Pop 
III star explodes in a 10^^ erg core-collapse supernova and we fol¬ 
low the non-uniform enrichment of the IGM and turbulent mixing 
that occurs as the blast-wave washes over neighboring halos. The 
simulation ends when metal-enriched gas is able to cool and col¬ 
lapse into dense, pre-stellar cores. The key results from this work 
are as follows: 

(i) The first metal-enriched stars form at z ~ 16.6 when a single 
Pop III supernova blast-wave washes over a neighboring mini-halo 
of comparable size that has yet to form its own Pop III star. Turbu¬ 
lence caused by the virialization of the halo and collision with the 
blast-wave allows metals to reach the interior of the halo. The gas 
within 0.01 pc of the dense core is uniformly mixed to a metallicity 
of~ 2 X 10“® Zq. 

(ii) Because the metallicity of the dense core is less than the gas- 
phase critical metallicity (Zcr ~ 10“^'^ Z©), the thermal evolution 
of the gas as it collapses is similar to the metal-free case until dust 
cooling becomes efficient at n ~ 10^^ cm“^. At this point, a sharp 
drop in temperature occurs that induces significant fragmentation 
into many small clumps with masses of roughly 10“^ M©. We run 
a second version of this simulation without dust and observe greatly 
reduced fragmentation. 

(iii) The velocity structure of the gas is dominated by turbulent 
motions at all scales, from ~ 100 pc (the outskirts of the halo) down 
to a few AU. The turbulent velocity, defined as the variance in the 
total velocity magnitude, roughly follows the local sound. At large 
radii, where efficient mixing occurs, the turbulent velocity reaches 
as high as Mach 2. The radial velocity also follows the local sound 
speed, but is somewhat smaller. The gas shows virtually no coher¬ 
ent rotation on any scale. 

(iv) By the end of the simulation, two Pop III supernovae have 
exploded in separate regions. The region in which the first super¬ 
nova explodes has relatively few halos. As a result, only the original 
host halo of the Pop III star exists as a potential site for second gen¬ 
eration star formation in this region. However, 15 halos exist in the 
region around the second supernova that have metal-enriched gas 
dense enough to form stars by the time that the larger structure en¬ 
compassing them has assembled at z = 10. Not counting the Pop III 
host halo, these halos have metallicities ranging from slightly over 
10“^ Z© to nearly 10“^ Z©. Thus, the first Pop II stars may have 
had a wide range of metallicities. 

We have identified one pathway, however small, from the ex¬ 
oticism of the early Universe to the normalcy of the present day. It 
is likely that multiple scenarios contributed to the transition from 
Pop III to Pop II star formation. While its overall contribution to 
Pop II star formation is probably minor, the external enrichment 
mechanism offers the best explanation for the origin of the most 
metal-poor stars observed and the clearest window into IMF of the 
first stars. 
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